Constraining the Equation of State of Dark Energy with Gamma Rays 
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Starlight in the Universe impedes the passage of high energy (e.g. TeV) gamma rays due to positron- 
electron pair production. The history of this stellar radiation field depends upon observations of star 
formation rate which themselves can only be interpreted in the context of a particular cosmology. 
For different equations of state of dark energy, the star formation rate data suggests a different 
density of stellar photons at a particular redshift and a different probability of arrival of gamma 
rays from distant sources. In this work we aim to show that this effect can be used to constrain the 
equation of state of dark energy. The current work is a proof of concept and we outline the steps 
that would have to be taken to place the method in a rigorous statistical framework which could 
then be combined with other more mature methods such as fitting supernova luminosity distances. 



I. INTRODUCTION 

Observations of type la supernovae 0, Q, galaxy sur- 
veys 0] and the CMB Q all seem to suggest that the 
expansion of the Universe is accelerating. This fact is 
usually explained by the presence of some component of 
energy density called dark energy which has an equation 
of state w < —1/3. The precise equation of state of this 
energy component is unknown and cannot be completely 
determined by the data (see e.g. (H). The density of 
the dark energy component may be a constant, it may 
be growing or it may be shrinking. Addition of the in- 
formation on angular distances provided by the baryon 
acoustic peak which can be seen in both the two point 
function of luminous red galaxies at z ~ 0.3 and in the 
two point function of CMB temperature observed at red- 
shift z ~ 1100 constrains the nature of dark energy quite 
well if the equation of state is a constant Q , but if the 
equation of state w is varying with cosmological epoch, 
the constraints on the nature of w(z) are less certain and 
exotic equations of states where the energy density is 
increasing are still compatible with the data (although 
difficult to motivate from theories of particle physics). 

In this work, we intend to constrain the equation of 
state of dark energy by calculating the effect of the differ- 
ent corresponding expansion histories upon the opacity 
of the Universe to high energy gamma rays. 

Stars are formed continuously throughout the history 
of the Universe. Some of the first stars (the ones with low 
masses) are still burning after order 10 10 years, whereas 
stars with masses larger than 10M Q last less than a few 
million years on the main sequence. The larger stars are 
more short lived, but much more luminous and also bluer. 
The continual life cycle of stars throughout the history of 
the Universe creates a background radiation in addition 
to the CMB. There are far fewer stellar photons than 
CMB photons but they have a higher energy: The initial 
photons from stars range from the near infrared to the 
ultra violet while a lot of the starlight fails to exit the 
dust in galaxies and is re-emitted in the infrared. This 
Extragalactic Background (star)Light or EBL therefore 
creates a spectrum of diffuse radiation ranging in fre- 
quency from the ultra-violet down to the border of the 



infrared and microwave bands at which point the CMB 
spectrum starts to dominate. 

CMB photons have energies of about 10 -4 eV so pho- 
tons with energies around 10 16 eV (i.e. 10 4 TeV) will be 
able to pair produce electron-positron pairs upon scatter- 
ing with the CMB Q. This sets the mean free path for 
such highly energetic photons to be very small in cosmo- 
logical terms, only around 10 kpc. We certainly therefore 
would never expect to see any such photons arriving from 
extragalactic sources and indeed photons with such high 
energies have never been positively identified (they could 
in principle form some fraction of ultra high energy cos- 
mic rays but the constraints on the photon fraction is 
rather tightly constrained The situation could be 
more interesting if axion like particles exist Q) 

Conversely, TeV energy gamma rays photons have been 
observed, and their arrival direction has been associ- 
ated with various extragalactic sources, including some at 
truly cosmological distances, (e.g. 3C 279 at a rcdhsift 
of z = 0.536) 0, [H. This is not completely surpris- 
ing since although TeV photons can instigate electron- 
positron pair production upon interacting with more en- 
ergetic starlight photons in the EBL, the number density 
of photons from stars is much lower than those in the 
CMB 0. By observing distortions in power law spec- 
tra from Blazars, several Cerenkov Telescopes have been 
able to place constraints upon the opacity of the Uni- 
verse [l(J [HI ■ This tallies well with observed star forma- 
tion rates and the standard cosmological model ACDM. 
However, if one assumes that the expansion history of 
the Universe is that corresponding to a varying equation 
of state for the dark energy component, the relationship 
between the observations of star formation made at dif- 
ferent redshifts and time is altered in such a way that the 
density of photons per square centimetre along the path 
of an incoming TeV gamma ray changes. Some equations 
of state make the Universe less opaque, some more. We 
will use the observed opacity of the Universe to place con- 
straints on expansion histories and show that this study 
can rule out some cosmologies which are otherwise con- 
sistent with the data. 

This extends our previous work looking at constraints 
upon Void universe cosmologies, cosmologies where the 
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apparent acceleration of the Universe is due to our be- 
ing located in a large underdensity [l2| . Such models are 
already very tightly constrained by observations, so our 
finding a new way of putting them under pressure using 
gamma rays is perhaps not of interest to everyone. Con- 
versely, dark energy with a varying equation of state is 
completely consistent with all available data sets, so any 
new constraints will be of interest and we hope to present 
one in this paper. 

First we will explain how we calculate the extra galac- 
tic background light from a given star formation rate, 
then we will describe how to calculate the opacity that 
light creates for high energy gamma rays. Then we will 
explain how we fit the star formation rate for different 
cosmologies before describing the procedure we use to 
obtain our results. After having presented our results we 
will become self critical and point out the various im- 
provements that need to be made to the different parts 
of the anaylsis in order for the method to become more 
rigorous and useful before concluding. The current work 
is therefore meant to be a proof of concept rather than 
providing a constraint which can be trusted at this stage 
as much as the standard fits. However we hope to con- 
vince the reader that the method outlined in this work 
could without too much work become a new tool for un- 
derstanding the nature of dark energy. 



II. THE EXTRAGALACTIC BACKGROUND 
LIGHT 



Because the extra-galactic background light is made 
of starlight, and we are surrounded by stars, it is quite 
difficult to see. Nevertheless, there are lower limits on 
its magnitude across the whole range of frequencies due 
to galaxy counts and positive detections from FIRAS in 
the infrared (for references, see the caption of figure [lj. 
It is possible to reconstruct the evolution and history of 
this radiation field using the star formation history of 
the Universe, which has been observed at various wave- 
lengths 0. 

We do this following closely the simplified approach of 
Finke et al. [Tij with a particular numerical implemen- 
tation which lends itself to different expansion histories. 
Let us outline our way of implementing those equations. 
The photon density of a blackbody is 



per unit time from a star of mass m and age t t 



n*(e, m, <*) 



8- 



dN 



dedV X 3 C exp [e/9(m, £*)]-! 



(1) 



where e = hu/m e c 2 is the dimensionlcss photon energy in 
units of electron mass. The dimensionlcss temperature is 
0(to, t*) = kBT{m,t Sf ) /m e c 2 where T(m,i*) is the tem- 
perature of a star of mass m and age i*. We are going 
to assume that this is single valued, i.e. we are going to 
neglect the effect of metallicity variation over time. We 
want to know how many photons of energy e are emitted 



N(e,m,U) 



dN 
dedt 



Ai:R{m, t*) 2 cn*(e, m, t r ) (2) 



where R(m, t*) is the radius of a star of mass m and age 
t» and c is the speed of light. We obtain the tempera- 
ture and radius R of a star of mass m and age i* we 
follow Finke et al. in using the results from the paper 
by Eggleton, Fitchett and Tout fla as well as the correc- 
tions outlined in the Finke work [LJ]. For a population 
of stars of age t* but with different masses, we need a 
weighted integral over the range of masses to include the 
effects of all the different stars. We therefore define the 
specific photon generation rate (per unit solar mass) 



N{e,U 



dN 
dedt 



£(m)N(e,m,t*)dm (3) 



where £ (m) is the initial mass function (IMF) of the stel- 
lar population. We will assume that this also doesn't vary 
over time. We assume a Salpeter initial mass function for 
the stars that are produced, namely £(m) oc to -2 35 with 
stars being produced across a range of m m i n = Q.IMq < 
in < m max = IOOMq. Here we have also introduced the 
effects of dimming, the parameter /(e) is the fraction of 
photons which do not make it out of galaxies from the 
blue end of the spectrum due to dust. The function /(e) 
that we have adopted is the one set out in [l4| although 
we have tried other dimming choices without changing 
the conclusions too much. 

The absorbed radiation is re-emitted in the infra red 
and the total integrated luminosity of this dust per solar 
mass of a stellar population of age t* is 



Ld{t*) = TUeC 2 



1 



V/W 



-1 )eN(e,t*)de (4) 



The energy emitted per solar mass of a particular stellar 
population created at time tf orm over its lifetime up to 
the time of the observer (e.g. today) to is 



J* (c, ^0 ) tform) — 

*° N[(l + z(t))e,t-t fo 

tform 



m e c 2 e 



(l + z(t)) 



l -dt (5) 



The dust re-radiates the absorbed energy thermally. 
There are at least three populations of different kinds 
of dust re-emitting the absorbed radiation at different 
temperatures so that we have 



jd(c, to, tf orm ) 



15/ n £d(i — tform) 

t form ti ^e* (n) (i + z(t)) 



E 



[(l + z(t))e 



cxp[(l + z(t))e/Q d{n) ] - 1 



dt (6) 



The three fractions and temperatures corresponding 
to the three different populations of dust are fitted by 
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©d(n) 


fn 


1 


7 x 10" 9 


0.7 


2 


25 x 10" 9 


0.15 


3 


76 x 1CT 9 


0.15 



TABLE I: Values for the fraction and temperature of the three 
populations of dust responsible for re-emitting the absorbed 
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FIG. 1: T/ie _B_BL spectrum produced by our code vs. the 
data for the ACDM case (solid green curve). This is plotted 
against actual data where the EBL light has been detected (blue 
circles) and other constraints which only place lower limits 
(green triangles). References for the data points are set out 
in the caption of figure 4 of Jla l. 



hand and the good fit values we obtain are listed in table 
HI The precise values of these numbers is not critical as 
it is in general the blue part of the spectrum which is 
responsible for the dimming of the gamma ray photons, in 
other words as far as opcity is concerned, the magnitude 
of the dimming of the bluest photons is more critical 
than the frequency at which they are re-emitted in the 
infrared. The overall spectrum of energy density at to is 
finally given by 



eu E BL(t,t ) 



; / p[z{t)]{j*(e,t ,t) +j d (e,t ,t)} dt 
>t star t 

(7) 

where p(z) is the star formation rate in units of solar 
masses per unit time per unit volume. To compare with 
observations, the energy density can be converted into 
intensity in units of nW m~ 2 sr -1 



el(e,t ) = —€u E BL{e 7 t ) 
Air 



(8) 



and the results of this spectrum can be seen for typ- 
ical parameters in figure [T] This therefore is the 
method by which the extragalactic background radiation 
is calaulated as a function of redshift. Next we need to 
see how this radiation field creates a finite mean free path 
for high energy gamma rays. 



III. SCATTERING OF HIGH ENERGY GAMMA 
RAYS WITH STARLIGHT 

To calculate the opacity of the Universe to photons we 
use the expression for the inverse mean free path (see e.g. 



8E' 2 

^ gamma J e 



n(e) 



sa(s)dsde 



(9) 

where s min = (2m e c 2 ) 2 , e min = s min /(AE y ) and er(s) is 
the cross section for electron-positron pair production 



tr(a) =mi\ (l-/3 2 )(3 - 4/3 4 )ln 



1 + /3 



1~(3 



2/3(2 - /3 2 



(10) 

where y = s/m^c 1 is the dimcnsionless centre of mass 
energy, P = y/l - 2/y and or is the Thomson cross sec- 
tion. 

We then calculate the opacity as a function of redshift 
by integrating 



dr —c 

dz A 77 (l + z)H(z) 



(11) 



then the attenuation in photons (i.e. the probability of 
them getting through from a given time t) is then given 
by P = exp(-r(i)). 



IV. FITTING THE STAR FORMATION RATE 

The star formation rate has been observed as a func- 
tion of redshift in the form of observations of galaxies at 
various redshifts [l3j]. In order to interpret these results 
in terms of a particular cosmology, one has to make a 
rescaling which reflects the different luminosity distances 
between the new cosmology and the reference cosmology 
in which the results are reported [l3| 



L(z) 



d 2 c (z) 



V c (z,Az) d 3 c {z + Az) - d 3 c (z - Az) 



(12) 



In this equation, the second term reflects the fact that 
the star formation rate is proprotional to the inverse of 
the comoving volume V c {z) corrcsponsding to the redshift 
bin between z + Az and z — Az while the second term is 
obtained taking into account the proportionality between 
Luminosity and comoving distance squared L cx d 2 c . 

Having re-scaled the star formation rate data and its 
errors bars we then fit it with the ansatz 



bz 



p* = 



1 



[z/c) d 



(13) 



and we fit the values of a, 6, c, d using a monte carlo 
method to find a best fit. This leads to some random 
fluctuations in the final results due to degeneracies in fits 
of a, &, c, d to the data however, as we shall see, the overall 



trend is clear. We have tried using the median value of 
the fit rather than the best fit but it doesn't significantly 
reduce the noise in the final data. Use of a larger star 
formation rate data set in future might help in reducing 
this scatter. We shall come back to this issue when we 
discuss the issues with errors that need to be addressed 
in the future in order for this to be a reliable method of 
constraining dark energy in section IVII 



V. EXACT PROCEDURE AND RESULTS 

This section describes the procedure for obtaining re- 
sults. First we choose a cosmology which has dark en- 
ergy with a nonstandard equation of state. The way we 
parametrise the equation of state is the following:- 



w(z) 



w 



1 



(14) 



We assume that the cosmology is flat and that £Im = 
0.276 which is close to the WMAP value Q. We also 
assume a Hubble parameter of h = 0.7. Ideally we would 
like to relax these assuptions and allow these parameters 
to vary, but the difficulty we will encounter in coming 
up with a reliable \ 2 statistic for the dimmming of pho- 
tons in different cosmologies makes this difficult. We will 
discuss these issues in section IVll 

Having chosen a particular cosmology (i.e. a particu- 
lar pair of w a and w a ), we obtain the comoving distance 
as a function of rcdshift as well as the time vs. red- 
shift relationship. We then apply the appropriate rescal- 
ing outlined in equation (Tj"2"j) to the data to make 
up for the different between this particular cosmology 
and the ACDM cosmology assumed in reference [l3| • We 
then vary the parameters a, 6, c and d in order to fit the 
rescaled star formation rate data using a monte carlo 
method, basically by using the Metropolis-Hastings algo- 
rithm and looking for the best fit value. Then we produce 
the extragalactic background light spectrum throughout 
the history of the new cosmology using the fitted values 
of a, b, c and d and the procedure described in section HI1 
Next we follow the procedure outlined in section Mil to 
work out how far high energy gamma rays can propagate 
through that newly derived radiation field. In particular 
we calculate the position of the line corresponding to the 
t = 1 opacity in the E-y — z plane (gamma ray energy vs. 
redshift). We then compare this line with the line which 
has been published in the literature from the Cercnkov 
telescopes. An example of two such lines, one of which 
is compatible with the observed dimming and the other 
which is not can be seen in figure [2] 

These r = 1 lines correspond to a one to one mapping 
between E 1 and z which we will denote E 1 (z)\ T=v We 
will call the reference limit on the opacity obtained by 
Magic f(z) = Ejmagic{ z )\ t= i while the function ob- 
tained for a particular w Q — w a cosmology we will call 
g{z,w ,w a ) = E iC osmology{z)\ t=v We can then de- 
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FIG. 2: Figure showing how some cosmologies are seemingly 
incompatible with observations. The black solid line is the 
constraint on the r = 1 gamma ray opacity from JTft l. The 
Green dot-dashed line and the red dashed line correspond to 
two different cosmologies with equations of state indicated in 
the figure (both are for h = 0.7 and Qm — 0.276). The cos- 
mology corresponding to the Green line is compatible with the 
gamma ray opacity data while that correspidning to the red 
line is not. 
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FIG. 3: The value of the overhang parameter $ defined in the 
text for different values of w and w a . It is clear that cos- 
mologies in the lower left part of the diagram are disfavoured. 



fine an overhang parameter f(m , w a ) as 



$>(w ,w a ) 



H(f(z)-g(z, w , w a )) (f(z) - g(z, w a , w a )) 

(15) 



where H{x) is the Heavysidc step function. This is sim- 
ply the square of the area that the r = 1 opacity curve 
overhangs the reference opacity curve on the E~—z plane. 
We plot this overhang parameter as a function of w and 
w a in figure [3l This figure shows clearly that cosmologies 
with low values of wq and w a are disfavoured by looking 
at the gamma ray opacity of such universes. This is in 
contrast to the cosmological data, which restricts the val- 
ues of w a and w a to lie inside a diagonal ellipse as can 
be seen in figure 0J Comparison of figures [3] and [5] shows 
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FIG. 4: Constraints on the parameters w and w a as obtained 
from Supernova Data, Baryon Acoustic Oscillations and Cos- 
mic Microwave Background Data in reference The two 
contours correspond to 0.99 and 0.95 times the maximum 
probability using supernova data from z = 0.03 upwards. This 
should be compared with figure [3| 



that regions which are not constrained by any cosmolog- 
ical measures are placed under pressure by the gamma 
ray opacity method outlined in this paper. The method 
therefore can therefore in principle regions of parameter 
space which the rest of the cosmological data set is con- 
sistent with. The current work is a proof of principle for 
this method. In the next section we will outline what 
needs to be done in order for this method to be placed 
on a stronger footing. 



VI. OBTAINING ERRORS AND FUTURE 
IMPROVEMENTS 

In the previous section we saw that there is tension be- 
tween the observed gamma ray opacity of the Universe 
and some combinations of the equation of state parame- 
ters w and w a which cannot be ruled out with the normal 
cosmological data. We also defined an overhang param- 
eter ^(w ,w a ) which measured how baldy each cosmol- 
ogy disagreed with the gamma ray opacity constraints. 
In principle it would be great to obtain a mapping be- 
tween '5 and a probability, then this probability could be 
brought together with the other coventional constraints 
on the dark energy equation of state in a combined anal- 
ysis. However this is not yet possible, for a number of 
reasons that we have until now not mentioned. In this 
section we will outline quite a number of these issues 
which will make up a to-do list which needs to be com- 
pleted in order for the method outlined in this work to 
be usefully employed. 

In order to turn the ^> statsictis into a x 2 statistic we 
would have to compare the overhang with errors on the 
gamma-ray opacity obtained at various rcdshifts. The er- 
ror on the constraint on the gamma-ray opacity provided 



by Magic is not quoted in their paper because they have 
not yet been able to calculate it in a robust way. They 
obtain the opacity as a function of rcdshift by looking 
at the distortions in apparently power law spectra which 
are interpreted as being due to electron-positron pair pro- 
duction as the photon crosses the intergalactic radiation 
field. However, this is an assumption which it is difficult 
to prove. Progress into coming up with a more rigorous 
statistic has been made in a recent paper by the Fermi 
collaboration working with more photons |18| ■ Hopefully 
the advent of the Cerenkov Telescope Array will allow 
the gamma ray observational community to come up with 
more robust meausurements of the gamma ray opcity as 
a function of redshift based upon more statistics and in- 
cluding errors. 

The method for calculating the extragalactic back- 
ground light is full of assumptions, for example, the shape 
of the initial mass function and its evolution over time. 
In our work we have assumed that one has a Salpeter imf 
throughout the entire history of star formation but this is 
clearly an over-simplification. Furthermore we have not 
taken into account the way that mctallicity changes over 
time and the variation of mctallicity in different spatial 
regions of the Universe (both of which of course have been 
observed to vary.) Only by measuring such variations and 
shapes as a function of redshift, or at least constraining 
them, then folding these variations and uncertainties into 
the analysis to understand the uncertainty in the extra- 
galactic background radiation field would we be able to 
obtain a reliable error on the dimming in a particular 
cosmology. 

Finally in the future we may be able to collect larger 
and more reliable data on the star formation histroy with 
smaller errors. This would enable us to increase our con- 
fidence in our fitting of the star formation rate parame- 
ters a, 6, c, d, indeed also to test this parametrisation and 
reduce a lot of the noise in figure [3] 



VII. SUMMARY AND CONCLUSIONS 

In this work we have suggested a new way of constrain- 
ing the equation of state of dark energy. We have calcu- 
lated the gamma ray opacity of Universes with different 
equations of state by rescaling the star formation rate his- 
tory according to each new cosmology. Some universes 
(i.e. pairs of values of w and w a in equation (|14p ) which 
are not ruled out by the normal cosmological data set 
would lead to an extragalactic background light in com- 
patible with the observed arrival of photons from distant 
objects. 

We have also outlined the steps which need to be taken 
to improve this technique and to place it on a firm statis- 
tical footing. This of course involves obtaining new and 
better data, both from gamma-ray telescopes constrain- 
ing the opciaty of the Universe better and from Galactic 
surveys learning more about the star formation rate of 
the Universe. 
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If these hurdles can be overcome it seems that the tech- 
nique presented in the current work could cut into regions 
of dark energy equation of state parameter space that 
cannot yet be probed be other means. 
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